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ABSTRACT 

We investigate the recent claim of ‘photon underproduction crisis’ by Kollmeier et 
al. (2014) which suggests that the known sources of ultra-violet (UV) radiation may 
not be sufficient to generate the inferred H I photoionization rate (Thi) in the low 
redshift inter-galactic medium. Using the updated QSO emissivities from the recent 
studies and our cosmological radiative transfer code developed to estimate the UV 
background, we show that the QSO contributions to Thi is higher by a factor ~2 as 
compared to the previous estimates. Using self-consistently computed combinations of 
star formation rate density and dust attenuation, we show that a typical UV escape 
fraction of 4% from star forming galaxies should be sufficient to explain the inferred 
Thi by Kollmeier et al. (2014). Interestingly, we find that the contribution from QSOs 
alone can explain the recently inferred Thi by Shull et al. (2015) which used the same 
observational data but different simulation. Therefore, we conclude that the crisis is 
not as severe as it was perceived before and there seems no need to look for alternate 
explanations such as low luminosity hidden QSOs or decaying dark matter particles. 
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1 INTRODUCTION 

Recently. iKollmeier et all d2014l . hereafter K14), used a cos¬ 
mological hydrodynamic simulation together with the lat¬ 
est measur ements of the H I colu mn density distribution, 
/(Am), bv iDanforth et all d2014h in the low -2 intergalac¬ 
tic medium (IGM) and reported a H I photoionization rate 
(Fhi) at 2 = 0. This is 5 times higher than the one (refer to 
as F^f 1 ) obtained from the theore tical estimates of c osmic 
ultraviolet background (UVB) by lHaardt fe Madaul J2012I . 
hereafter HM12). This apparent discrepancy has led to the 
claim of a ‘photon underproduction crisis’ suggesting that 
the origin of more than 80% of H 1 ionizing photons is un¬ 
known and perhaps generated from non-standard sources. 

For a given sight-line in a cosmological simulation, the 
inferred /(Am) depends on the assumed Thi, the distribu¬ 
tion of gas temperature and the clumping factor of the region 
producing the Ly-a absorption. The latter two quantities de¬ 
pend not only on the assumed initial power spectrum but 
also on various feedback processes that inject energy and mo¬ 
mentum into the IGM from star forming galaxies. Therefore, 
the Thi estimates using the /(Ahi) will depend on how re¬ 
alistic the various feedbacks used in the simulation are. K14 
hav e used the smoo th particle hydrodynamics code GADGET 
2.0 dSpringelll2005l ) that includes feedback from galaxies in 
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the fo rm of mo mentum driven win ds (lOppenheimer fc Pavel 
120081) . However. I Dave et all d20ld ) suggested that these feed¬ 
backs produce negligible effect on /(Ahi) for Am < 10 14 
cm -2 . 

Recently, IShull et al.l (120151 ) have independently esti¬ 
mated T hi, using the same observed data but simulated spec¬ 
tra obtained using the grid based Eul erian N-body hydrody¬ 
namics code enzo dBrvan et al.ll2014l ). They found a smaller 
r hi than K14 but it is still a factor 2 higher than T(| j vl . They 
attributed the decrease in the derived T hi as compared to 
K14 to the differences in the implement ations of fe e dback 
processes in the simulations used. While fshull et akl d2015l ) 
reduced the apparent tension, it still requires an appreciable 
contribution to the UVB from galaxies when one uses the 
previously estimated QSO emissivity. 

I11 this study, we revisit the UVB calculations at 2 ~ 0 
using the numerical radia tive transfe r code developed by us 
( Khaire fe Srianandl2013l ) in line with lFaucher-Giauere et aid 
( 20091 ) and HM12. We update the QSO and galaxy emissivity 
an d show that the Q SOs alone can provide the T hi inferred 
by lShull et aid d2015l ) and to get the Fhi inferred by K14, we 
need only 4% of the ionizing photons to escape from galaxies 
(and not 15% as suggested by K14) . Throughout this paper 
we use a cosmology with S1 a = 0.7, fl m = 0.3 and H 0 = 70 
km s _1 Mpc -1 . 
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2 THE RADIATIVE TRANSFER 

Following _ the _ stan dard _ proc edure 

(iMiralda-Escude fc Ostriked Il99fll: IShapiro et alJ 1 19941: 

Haardt fc MadauT 19961 : iFardal et al.l Il998[ Shull et all 

1999lb the average specific intensity, J vo (in units of erg 
cm" 2 s" 1 Hz' 1 sr' 1 ), of the UVB at a frequency vq and 
redshift zo is given by, 

1 r°° rll 

J M = ^ ^ dz (1 + zo ) 3 e*(*) e—(!) 

Here, (p is the Friedmann-Lemaitre-Robertson-Walker line 
element, e„(z) is the comoving specific emissivity of the 
sources and r e g is an average effective optical depth encoun¬ 
tered by photons of frequency uo at a redshift zo which were 
emitted from a redshift z > zo with a frequency v > i/q. The 
frequency v and vq are related by v = uo(l + z)/( 1 + zo). 
Assuming that the IGM clouds of neutral hydrogen column 
density, IVhi, are Poisson- distributed along th e line of sight, 
r e ff can be written as (see IParesce et aljfl980h . 

PZ POO 

Tes{yo,zo,z) m / dz' dA H i/(A H i, z)(l - e - ^'). ( 2 ) 
J zq Jo 

Here, /(JVhi,z) is the number of H 1 clouds per unit red¬ 
shift and column density interval having column density 
Ahi. The continuum optical depth t„/ is given by t v i = 
Nmcrm(i'') + AHeiCHei^') + A H eiio- H eii(^'). where, A; and 
< 7 i are the column density and photoionization cross-section, 
respectively, for species i and 1 / = ^ 0(1 + z')/(l + zo). 

We use the same /(Ahi,z) used by HM12 and neglect 
the contribution of He 1 to t v because of its negligible abun¬ 
dance at z < 6. We calculate r e fi following the prescription 
given in HM12. In the following section we provide the up¬ 
dated source emissivity. 


3 EMISSIVITY OF RADIATING SOURCES 

We calculate the UVB assuming only QSOs and galaxies are 
sources of the UV radiation. Therefore, e„(z) = e^(z) + e^(z) 
where (z) and e(?(z) are the comoving specific emissivity 
from QSOs and galaxies, respectively. 


3.1 Comoving QSO emissivity 

The e(?(z), in units of erg s _1 Hz -1 Mpc -3 , using the ob¬ 
served QSO luminosity function (QLF) at a frequency v is 
given by 

/ OO 

L v {z)tj>{L v ,z)dL v , (3) 

where, (^(L^ , z) is the QLF at z given in terms of specific 
luminosity L„ by 

= frl/L*) [ ( L v /L T 71 + (W^T 72 r\ (4) 

and using the absolute AB magnitudes, M, by 

4 >{ M ) = <^m[ 10°' 4(71+1)(m_m * ) + 10 °- 4 (t' 2+ 1 )( m - m, ')]- 1 ; 

( 5 ) 

where (j>*M = 0.921 We use L]= 0.01L* to calculate 
the ej?. Note that, for most of the QLFs at z < 3.5 the faint 
end slope 71 Js —1.4 (see Table 2), for which the value of 
is insensitive to the values of < 0.01L*. For example, 


the maximum difference in the calculated for L™*” = 0 
and 0.01L* is less than 5% for 71 > —1.4 and 72 = —3.5. 

In Fig.l ( left panel), we plot the L„(/>(Z/„) estimates 
against g-band magnitudes from various studies at two dif¬ 
ferent z. The area under each curve is proportional to 
the respective emissivity at g-band (e^). It is clear from 
the Fig.l that, as com pared to old QLF m easurements of 
iBovle et al.l (120001) and Croom et al.| I2004h . using the new 
measurements of Croom et al.| d2009l . hereafter C09) and 
IPalanaue-Delabrouille et al.1 ( 201 3 . hereafter PD13) will give 
a larger e®. This is indeed the case, as demonstrated in the 
right panel of Fig.l where we plot the for these QLF 
measurements. We have also plotted the tf converted from 
the ej 12 (z:) given at the H 1 Lyman limit (i.e at 912A) by 
HM12 using the relation log(e^) = log(e® 12 )-|-0.487 which is 
consistent with the spectral energy distribution (SEP ) used 
by HM12. This ej? i s consistent with IBovle et al.l (l2000l l and 
ICroom et al.l (1200414 which is smaller by factor ~ 1.5 to 2 as 
compared to C09 and PD13. 

In our study, we use the latest QLF measurements as 
summarized in Tabic [T] The first and second columns give 
the reference and the wavelength (Ab an d) at which the QLF 
is reported, respectively. For 0.3 < z < 3.5, in each redshift 
bin we fit the observed QLF with the form given in Eq. [5] 
using an idl MPFIT routine by fixing the values of 71 and 
72 to those reported in the respective references (our fits 
are also presented in Fig [T]). We use our best fit <\>* M and 
M* to obtain eAband(z) (see Table [TJ. At other redshifts, 
we take the best fit QLF parameters given in the respective 
references and calculate the e^ hand {z). In Fig 1 ( right panel), 
we show that the ef (z) at z < 2 obtained using our fit is 
consistent with the pure luminosity evolution (PLE) models 
of PD13 and C09. 


We convert the e^bandi 2 ) i n t° e gi 2 ( 2 ) us i n 8 the broken 
power law QSO SED L„ oc v~ a , which we adopt for our 
UVB calculations. In the soft X-ray re gime above energy 
0.5 keV (A ^ 2 4. 8 A) we use a = 0.9 (INandra fe Pourida 
Il994lh Following l&tevans et al.l d2014lb we use a = 1.4 for 
24.8 < A < lOOOAQand a = 0.8 for 1000 < A < 2000A. For 
A > 2000 A we use a = 0.5. 

For the observed QLF at z < 3.5, the SEDs used to 
perform continuum A-corrections in the original references 
(L„ oc v~ a ; a' is given in the last column of Table |T]) are 
different from our adopted SED at A ^ 1500A. For con¬ 
sistency, we recompute the specific emissivity £Arest( 2 ) a t 
Arest = Aband /(1 + z) using a' from the corresponding refer¬ 
ence and then use our adopted SED to convert tArest ( z ) ' rl to 
e 9 i 2 (z). This is not needed for z > 4 where the QLFs are ob¬ 
tained at 1450A in the QSO’s rest frame using appropriately 
matched filters without applying additional A-corrections. 
The errors on the ej 12 (z:) given in the Table 1 are the maxi¬ 
mum and minimum difference we get using the errors in 71 
and 72 given in original references except at z = 0.15. In 
this case the error on e ® 12 is the difference we get in ef 12 if 
we use = 0.1 and 0.001. 

All our e ^ 12 (z) measurements as a function of z are plot- 


1 IStevans et al. | d2014T) found a = 1.41 ± 0.15 using QSO com¬ 
posite spectrum down to about 500 A. We extrapolate it upto 
A f>j 25A. Note that the UVB at A < 500A has negligible contri¬ 
bution to Thi- 
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Figure 1. The L u <j){L u ) is plotted against M at g-band for two redshifts using various QLF reported in the literature. The area under 
each curve is proportional to the e® (left panel). The eg obtained for these QLFs, our eg from the fit to the PD13 and C09 QLF (cyan 
p oints; see tab l e 1 for details) and the e® ( z) inferred from e^ 0 (z) of HM12 are al s o plot ted (right panel). Here, the best fit PLE models 
of iBovle et al.1 (|2000h . fCroom e t al.1 (12 004). with 2SLAQ data of I Richards et ahl (|2005l h C09 and PD13 and the luminosity dependent 
density evolution model of lBongiorno et all (120071) are used. 


ted in Fig. [2] We fit these points using a functional form 
similar to that of HM12 and obtain the following best fit, 


4 2 ( z ) = 10 24 ' 6 (1 + z ) 5 


exp(— 0 . 362 ) 
exp(2.2z) + 25.1 


( 6 ) 


For comparison, in Fig(2j we show this best fit £gi 2 ( 2 ) along 
with the Cg 12 (z) used by HM12. For 2 < 3.5, our eg 12 (z) is 
higher than that of HM12 and the maximum difference of 
factor 2.1 occurs at 2 ~ 1.5. The peak in £ 912 ( 2 ) also changes 
from 2 = 2.2 from HM12 to 2 = 1.95. In Fig. [5] we also show 
the £912 ( 2 ) obtained using the PLE models of C09 and PD 13 
and the l uminosity e v olutio n and density evolution (LEDE) 
model of iRoss et all (120131 1. These are consistent with our 
fit in Eq|S] Note that the PLE models of C09 and PD13 
give identical values of ef(z) (see right panel of Fig.l) but 
differ slightly in the 4 2 ( 2 ) since they use different SED for 
continuum A'-correction (see Table 1). 


3.2 Comoving galaxy emissivity 

In iKhaire fe Srianandl (120141 . hereafter KS14), by match¬ 
ing the observed galaxy emissivity from multi-band, multi¬ 
epoch galaxy luminosity functions, we have determined self- 
consistent combinations of the star formation rate density 
(SFRD) and dust attenuation magnitude in the FUV band 
(Afuv) for five well known extinction curves. It has been 


found that the SFRD ( 2 ) and Afuv( 2 ) estimated using the 
average extinction curve of the Large Magellanic Cloud Su¬ 
pershell (LMC2) is consistent with various observations. 

Here, as our fiducial model, we use the £^( 2 ) computed 
from the SFRD ( 2 ) and Afuv( 2 ) obtained in KS14 for the 
LMC2 extinction curve (see tables 2 and 4 in KS14). Our 
SFRD at 2 < 0.5 is a factor ~ 3 higher than that of HM12 
(see also lMadau fe Dickinsonll2014l l. However, the difference 
decreases at higher 2 , and becomes less than 10% at 2 ~ 3. 
A small fraction, f eBC , of the generated H 1 ionizing photons 
(A < 912A) from the stellar population are assumed to es¬ 
cape through holes in galaxies (i.e by assuming that dust 
does not modify the SED at A < 912A). We assume that 
there are no He 11 ionizing photons (A ^ 228A) escaping the 
galaxy. This is a reasonable assumption in the 2 -range of our 
interest. We approximate the galaxy emissivity at A < 912A 
with a power-law e® oc i/ -1 ' 8 . The exponent is fixed to repro¬ 
duce the Thi obtained from the model spectrum itself. Note 
that the exponent and the total H 1 ionizing photons gen¬ 
erated inside the galaxy depends on the metallicity, initial 
mass function (IMF), stellar rotation rates a nd adopted evo¬ 
lutionary tracks (see iTopping fe Shull 2015lb In our galaxy 
models obtained from STARBURST99 llLeitherer et aljfl999l ). 
we use the Salpeter IMF with 0.4 times solar metallicity. 
See KS14 for a discussion on the uncertainties in estimating 
SFRD ( 2 ) and Afuv( 2 ) arising from the assumed metallicity 
and IMF. 
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Figure 2. The best fit comoving e^ 12 (z) used in this paper (blue curve) along with the c^ 12 {z) used by HM12. The green circles are the 
e ?i 2 gi ven in Table 1. The blue curve is simply a f it to these points. The €q 12 (z) obtained using the PLE model of C09 (dot dash curve) 
and PD 13 (dotted curve) and the LEDE model of I Ross et all (120131 , orange curve) is shown. 


In addition to this, we have also included some of the 
diffuse emission from the IGM clouds. We model the He II 
Ly-cr and He II Balmer continuum recombination emission 
following the prescription given in HM12 and the Lyman con¬ 
tinuum emission due to reco mbination o f H I and He I I using 
the approximations given in IFaucher-Giguere et al.l (120091 '). 
We do not include the contributions to UVB from He I re¬ 
combinations and the two photon continuum. These contri¬ 
butions are negligibl e, and i f included, can increa se Thi by 
a maximum of 10% (IFaucher-Giguere et al.l 120091 ') . We also 
do not include the resonance absorption of Hell which has a 
negligible effect on Thi, especially at low -2 (see HM12). 


4 RESULTS AND DISCUSSIONS 


Here, we focus on the H 1 photoionization rate, Fhi, obtained 
using our UVB model. This is defined as 


Fhi = 



du 


47 rj v 
hv 


CHl(i') , 


( 7 ) 


where I'm corresponds to A = 912A. In Fig. [3] we summa¬ 
rize various available Fhi measurements as a function of 2 . 
In particular, denoting Thi,13 = Thi x 10 13 s~ , the points 
of interest for the present study are Fhi,13 ~ 1.8 at 2 = 0.1 
as inferred bvlKoll mcier et ah ( 2014 ) which led to the claim 
of a ‘photon underproduction crisis’ and the very recent es¬ 
timates of Fhi, 13 (g) = 0.46(1 + z) 4 ' 4 at 2 < 0.5 found by 
I Shull et al] (1 20151 ). 


To begin with, we validate our code by reproducing the 
results of HM12. In Fig [3] we plot the Fhi(;z) determined 
by HM12 (long dashed curve) and the result of our code ob¬ 
tained using the e®( 2 ), SED and f esc used by HM12 (dotted 
curve). Both match with each other within ~ 5% accuracy. 
The minor differences noticed can be attributed to the dif¬ 
ferent metallicities used and contributions of some of the 
diffuse emission processes ignored in our model. Having val¬ 
idated our code, we use the updated QSO emissivity £ 912 ( 2 ) 


(see Eq. H and the ( 2 ) mentioned above (in Section 3.2) 
to calculate the UVB (and hence Thi) for different values of 

/esc- 


When we use only the QSOs as the source of the UVB 
(by taking f esc = 0) and use our updated e„(z), we get the 
Thi at 2 < 0.5 very c l ose (i .e wi thin 10%) to the v alues pre¬ 
dicted bv IShull et al] J2015I ) and lShull et al.l (120141). We find 
that t he one-sided ionizing flux 4>o, as defined in Shull et al.l 
(120151 ) is to be 5030 cm -2 s _1 for our UVB at 2 = 0 as com¬ 
pared to 5700 cm -2 s _1 obtained bv IShull et al.l d2015l ). How¬ 
ever, because of the statistical uncert ainties in the obse rved 
/(Am), the Thi and 4>o predicted by 1 Shull et al.l (120151 ) can 
be even higher. Our Thi, 13 values are 0.41, 0.94, 1.9 and 3.3 
at 2 = 0, 0.2, 0.4 and 0.6, respectively. These are ^2 times 
higher than the corresponding T ([( I values. Now, instead of 
using our e„{z) fitting form, if we take the best fit PLE 
models given in C09 and PD13 (see Fig[2j for 2 < 2.2, and 
estimate the UVB by assuming eg 12 ( z ) = 0 at 2 > 2.2, we get 
the Thi, 13 at 2 = 0 to be 0.48 and 0.39, respectively. It shows 
that, irrespective of our QLF fits and the fitting form, the 
updated QSO emissivity will lead to Thi ~ 1.7 t o 2.1 x r™. 
There fore, we conclud e that the Thi inferred bv IShull et al.l 
d2014l ) and I Shull et al.l (120151 ) can be explained by the QSOs 
alone without requiring any significant contribution from the 
galaxies (i.e with f esc = 0). This is consistent with many low- 
2 u pper limits on average f esc measured in samples of galax¬ 
ies (ISiana et al.l 20071 : ICowie et all 120091 : ISiana et al.ll2Q10l : 
iBridee et al.l 20ld : Leitet et al.ll2013l ). Therefore, there is no 
real photon unde rproduction c risis w hen we consider the Fhi 
measurements of lShull et al. I 1201 A) . 


In our UVB calculations with / eac = 0, we use a different 
QSO SED and an updated e?(g) as compared to HM12. 
However, since the Fhi oc (3 + a) -1 , changing a from 1.57 
(HM 12 ) to 1.4 at A < 912A increases the Thi by only 4%. 
The main difference in Thi between our UVB and that of 
HM12 arises because of the updated e„(z). It is important 
to realize that even though the e ® 12 used by us matches with 
e 9 i 2 °f HM12 at 2 = 0, the local UVB is contributed more 
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Figure 3. The Thi vs z obtained for our UVB with / esc = 0 ( solid curve) and with f esc = 4% ( dot-dash curve) along with the Thi from 
HM12 ( dash curve) is plotted. The dotted curve shows t he Tht when we obtain the U VB using our code with the ,, (z ) . SEP and f esc 
taken f rom HM12. The Thi me asurements at high-z bv lFaucher-Giguere et ahl (12008! ) ( squares). bvlBo lton &; Haeh nclt ( 20(53) ( circles) 
and bv lBecker fc Boltonl 120 I 'A) ( triangles) are shown. At low-z, the lower limit fm bvlAdams et al.l fl20 1 ll ) using non-detection of Ha 
from UGC 7321 ( arrow, for more details on the validity of it see sec. 3.2 of [Shul]_eGalJ|201^), the Thi whi ch is found consistent w i th the 
cosmic metal abundances bv IShull et al.l (1 201 fl i ( diamond ) and the inferred Tm of IKollmeier et al. ( 2014 1 (star) and lShull et al.l (1201 Oi l 
( green curve with diamonds) are also plotted. 


by ionizing photons coming from high-z, up to z ~ 2 , where 
the mean free path for H I ionizing photons is very large and 
6 ^ 12 ( 2 ) peaks. 

Next we explore the f esc requirements i n orde r to repro¬ 
duce the Thi inferred by iKollmeier et al.l (120141') . For sim¬ 
plicity we run models keeping / esc constant over the full z 
range. We find / es c = 4% is needed to get the Thi, 13 = 1.8 
at z = 0.1 (see Fig [3]). Interestingly, the / esc needed in our 
calculations is much less than the / esc = 15% required in the 
HM12 UVB model. Apart from 2 times higher QSO emissiv- 
ity, it is partly because of our ~ 3 times higher low-z SFRD 
as compared to HM12. Note that, the value of / eac ~ 0.02% 
used by HM12 to estimate UVB at z = 0.1 is extremely 
small compared to various low-z observations. In passing, 
we note that for our models with different combinations of 
SFRD and Afuv explored for different extinction curves in 
KS14, we require / esc values similar to or less than what we 
have obtained here for our fiducial model. 


In order to compare with observations, we use a relative 
escape fraction, / es c ,rei, defined as / esc ,rei = /esc x 10 O 4Afuv . 
To match the Thi of K14, the model of HM12 that assumes 
Afuv = 1 at z < 2, will require f es c , ie \ = 38% while we need 
only /esc.rei = 15% for our fiducial LMC2 model at z = 0 
(where we determined Afuv = 1.42). This / es c,rei = 15% is 
about a factor ~ 2 higher than the low z upper limits on 
/esc.rei given in various studies of galaxy samples as men- 
tion ed above. However the f eac observed in individual galax¬ 
ies ( |Borthakur_et ah 201411 and many t heoretical estimates 
('e.g. lKimm fe Cenll2014l : lRov et al.ll2014l l are consistent with 
it. Therefore, we conclude that with the updated QSO and 
galaxy emissivities presented here, even if we wish to gener¬ 
ate Thi inferred by Klf, the required f eS c of ionizing photons 
from star forming galaxies is not abnormally high enough to 
warrant an alternate non-standard source of the UVB. 


Interestingly, our updated QSO emissivity alone 
can reproduce the Th i measurements at high z 
dBolton fc Haehneitl 120071 : iBecker fc Boltonl l2013h up to 
z ~ 2.7. However, f esc = 4% gives a Fhi(2 ) which marginally 
overestimates the Thi measurements at 2 < z < 3 (see 
Fig.3). Irrespective of the low-z Thi, at z > 3 one needs 
galaxies to contribu te more to the UVB (however see, 
iGiallongo et al.ll2015l l. At higli-z, using the observations of 
H I and He II Ly-a forest, it will be possible to const raints 
the /esc from galaxies (see, iKhaire fc Srianandl l2013l f. We 
plan to do this in the near future. 


5 SUMMARY 

The recent claim of a ‘photon underproduction crisis’ 
(IKollmeier et al.ll2014lf requires the low-z Thi to be 5 times 
higher than the one obtained by the UV B mode l of HM 12. A 
similar investigation performed by I Shull et al.l (120151 1 finds 
a lower Thi which is still 2 times higher than that of HM12. 
Here, we present an updated H I ionizing QSO emissivity 
by using recent QLF measurements. It turns out that this 
emissivity is a factor of 1.5 to 2 times higher than what is 
used by HM12 at 0.5 < z < 2.5. We estimate the UVB using 
this emissivity with the help of a radiative transfer code de¬ 
veloped by us . We show that Q SOs alone can give a factor 
2 required bv IShull et al.l (12015 1. Using our updated SFR.D 
which is ~3 times hig her than HM12 at low-z, to get the 
Fhi predicted bv lKollmeier et ahl (120141 ') we require only 4 % 
of the ionizing photons generated by galaxies to escape into 
the IGM. Therefore, there is no need to look for additional 
sources of ionizing photons such as hidden QSOs or decaying 
dark matter particles. 
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Table 1 . Details of observed QLF used to get e^ 12 in our study. 


Reference 

•^band 

z 

logtp* 

M* 

71 

72 

lo § e A band 

loge^ 2 

a' 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

Schulze et al. (2009) 

4450A 

0.15 

-4.81 

-19.46 

-2.0 

-2.82 

24.30 

23.83±0.20 

0.5* 

Croom et al. (2009) 

4686A 

0.54 

-5.98 

-24.10 

-1.4±0.1 

-3.5±0.1 

24.59 

24.14±0.05 

0.3 

PD13 f 

4686A 

0.87 

-5.74 

-24.68 

-1.4±0.1 

-3.5=t0.1 

25.07 

24.59±0.05 

0.5 



1.25 

-5.81 

-25.65 



25.38 

24.90±0.05 




1.63 

-5.80 

-25.80 



25.45 

24.98±0.05 




2.01 

-6.00 

-26.71 



25.62 

25.17±0.05 


Ross et al. (2013) 

7480A 

2.25 

-5.83 

-26.45 

i q + 0.5 
± ‘°- 0.2 

q c + 0.3 

25.64 

25 . 06 t°-ii 

0.5 



2.35 

-5.98 

-26.55 



25.53 

24 95”*” 0 ' 11 
z^.yo_ 0 re 




2.45 

-6.12 

-26.81 



25.50 

24 91 + 0 - 11 
0.16 




2.55 

-6.12 

-26.77 



25.47 

24.89+o'ig 




2.65 

-6.21 

-27.06 



25.51 

24.921°;“ 




2.75 

-6.26 

-27.27 



25.54 

24.96lg;“ 




2.90 

-6.29 

-27.22 



25.49 

24 Q 1+ 0 * 11 
^ , yi - 0.16 




3.12 

-6.38 

-27.13 



25.36 

24 79+ 011 
^•‘ y - 0.16 




3.37 

-6.66 

-27.63 



25.29 

24 72”*” 0 ' 11 

1 - 0.16 


Glikman et al. (2011) 

1450A 

4.0 

-5.89 

-24.10 

-1 6 -1 " 08 

-3.3±0.2 

24.80 

24 62”*" 0,50 
z^t.Dz_o so 

NAi 

Masters et al. (2012) 

1450A 

4.25 

-7.12 

-25.64 

-1.72±0.28 

-2.6±0.63 

24.32 

24 13+ 0 - 50 
^• iO -0.22 

NA 

McGreer et al. (2013) 

1450A 

5.0 

-8.47 

-27.21 

-2 03+ 0 - 15 
z - uo - 0.14 

-4.0 

23.78 

23.60l°“ 

NA 

Kashikawa et al. (2015) 

1450A 

6.0 

-8.92 

-26.91 

-1 92 -1-0 ' 24 

- L-y - 0.19 

-2.81 

23.13 

22.94±0.15 

NA 


Column (4) gives <j >* hi units Mpc —3 mag -1 , column (5) gives M* in AB magnitudes and column (8) and (9) gives ej? in units erg s —1 
Hz -1 Mpc -3 . The Aband at 1450 A, 4450A, 4686A and 7480A correspond s to FUV, B, g and i band, respecti vely. *We assume a' = 0.5 
consistent with the k -correction of I Schulze et alj (120091 ). ^PD13 stands for lPalanque-Delabrouille et al.l d2013l) . ^NA indicates that the 
K-correction is not applied. 
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